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Abstract The structure of compact star core filled by gluon matter plasma 
CI is investigated within the fluid-like QCD framework. The energy momentum 

tensor, density and pressure relevant for gluonic plasma having the nature 
of fluid bulk of gluon sea are derived within the model. It is shown that the 
©h| model provides a new equation of state for perfect fluid with only a single 

parameter of fluid distribution 4>(x). The results are applied to construct the 
equation of state describing the gluonic plasma dominated compact star core. 
i-C ■ The equations of pressure and density distribution are solved analytically for 

small compact star core radius. The phase transition of plasma near the core 
surface is also discussed. 

> 

' Keywords quark-gluon-plasma ■ fluid QCD ■ compact star • equation of 

state 

' 1 Introduction 

Recent experiments in the last decades on relativistic nuclear collisions shed 
light on the phenomena of hot plasma formed by dense quarks and gluons. 
Those experiments suggest that the quark gluon matter behaves more like a 
^> \ deconfined quark-gluon plasma (QGP) liquid [1-3]. 

This fact immediately encourages some models based on either the (rela- 
tivistic) hydrodynamic [4,5], or pure Quantum Chromodynamics (QCD) ap- 
proaches [6]. Within the hydrodynamic framework, the plasma is dominated 
by either quark [5] or gluon matter [4]. In particular, in the quark matter 
dominated plasma dissipative ideal hydrodynamics has been used to fit some 
experimental data at high energy heavy ion program at the Relativistic Heavy 
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Ion Collider (RHIC) [7] . The successful fit requires the models to take into ac- 
count very small value of the ratio shear viscosity over entropy [8-13]. However 
the puzzle must still be confirmed by the next coming experiments at the Large 
Hadron Collider (LHC) [14]. While the gluon matter plasma is motivated by 
the discoveries of jet-quenching in heavy-ion-collision at RHIC indicating the 
shock waves in form of March cones [15, 16]. 

On the other hand, in pure QCD, QGP is described as a quark soup before 
hadronization which is a phase of QCD, and exists at extremely high tem- 
perature and/or density. It is argued that this phase consists of almost free 
quarks and gluons. Therefore, the phase transition from the deconfined QGP 
to the hadronic matters or vice versa gets particular interest in this approach. 
It unfortunately turns out to the many body problems with large color charge 
which cannot be calculated analytically using perturbation. As a result, the 
main theoretical tools to explore QGP within QCD is lattice gauge theory. The 
lattice calculation predicts that the phase transition occurs at approximately 
175 MeV [17,18]. Many works on both analytical and numerical calculations 
have been done at a sophisticated level providing some interesting predictions 
like color glass superconductivity (CGS), color- flavor locked (CFL) phase and 
so on. 

Anyway, describing the QGP as either quark or gluon matter dominated 
plasma within hydrodynamic or pure QCD should be understood as conjec- 
tures in the world of QGP. The important point is, since the QGP contains 
many quark-anti-quarks and gluons, it is considerable to treat it using the 
well-established QCD. At the same time from hydrodynamic point of view, 
the experiments strongly suggest that QGP behaves like a fluid. In particular, 
in the scenario of viscous gluonic plasma this is required to form and propa- 
gate shock waves [19-21]. Moreover it dissolves into an almost perfect dense 
fluid [22]. Therefore it is plausible to describe it as a strongly interacting fluid 
system. In this sense, there are approaches based on unifying or hybridizing 
the charge field with flow field [23-28] . Recently, some works have constructed 
the models in a lagrangian with certain non-Abelian gauge symmetry to the 
matter inside the fluid [29,30]. 

In this paper, we follow the scenario of gluon matter plasma as suggested 
by hydrodynamic approach, while it is still governed by QCD. In this frame- 
work, the strongly interacting system of QCD is considered as a macroscopic 
fluid system rather than a result of many particles interactions. It is shown 
that one can derive the total energy momentum tensor of gluon plasma within 
the model. Further, the equation of state (EoS) relevant for plasma dominated 
compact star core interiors is constructed. In the model, the density (p) and 
pressure (p) are determined and related each other through the fluid field dis- 
tribution <f>(x). Hence, the EoS can be solved through the ordinary differential 
equation (ODE) of p to show the phase transition within the model. 

The paper is organized as follows. First we briefly introduce the underly- 
ing model of gauge invariant fluid lagrangian and discuss the relevant physical 
scale and region within the model. Then, the energy momentum tensor, den- 
sity and pressure in the model are derived and investigated. Subsequently it 
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is followed by constructing the relevant EoS for compact star core interiors. 
Before summarizing the paper, the analytical solution for the pressure and 
density distributions of small compact star core radius is given and discussed. 



2 The model 

Let us adopt the model developed by Sulaiman et.al. [30,31]. The model pro- 
poses to describe the QGP as a strongly interacting gluon sea with quarks and 
anti-quarks inside. The model deploys the conventional QCD lagrangian with 
SU(3) color gauge symmetry, that is, 

C = i&fdpQ - m Q QQ - \s«„S a » v + g s j;U a » . (1) 

Here Q and {/ M represent the quark (color) triplet and gauge vector field. g s 
is the strong coupling constant, J° = QT a j^Q and T a 's belong to the SU(3) 
Gell-Mann matrices. The strength tensor is S£ v = d^U? - d v U% + g s f abc U b U^ 
with f abc is the structure constant of SU(3) group respectively. It should be 
noted that the quarks and anti-quarks feel the electromagnetic force due to 
the U(l) field A^, but the size is suppressed by a factor of e/g s = y/a/a s ~ 

oiio- 1 ). 

Following the original model [30] , the gluon fluid is put to have a particular 
form in term of relativistic velocity as, 

E£ = (Etf,U o ) = «-0, (2) 

with m° = 7 v o(l,v°) and j v a = (1 — |v a | 2 ) . <j> = 4>{x) is a dimension one 
scalar field to keep correct dimension and should represent the field distribu- 
tion. It is argued that taking this form leads to the equation of motion (EOM) 
for a single gluon field as follow [30] , 

^ (7v«VV) + V ( 7v „0) = -g s j> dx (J a + F£) , (3) 

where J* is the covariant current of gluon field, and F£ is an auxiliary function 
which can be found in the original paper [30] . It has been concluded that Eq. 
(3) should be a general relativistic fluid equation, since at the non-relativistic 
limit Eq. (3) coincides to the classical Euler equation. 

More precisely, this fact provides a clue that a single gluonic field 
may behave as a fluid at certain scale, beside its conventional point particle 
properties with a polarization vector e M in the form of U£ = <j>. One can 
then consider that there is a kind of "phase transition" , 

hadronic state < — > QGP state . (4) 
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As the gluon field behaves as a point particle, it is in a stable hadronic state 
and is characterized by its polarization vector. On the other hand in the pre- 
hadronic state (before hadronization) like hot QGP, the gluon field behaves 
as a highly energized flow particle and the properties are dominated by its 
relativistic velocity. 

One should also recall that the wave function U£ for a free massless particle 
satisfies d 2 U^ = with a solution U" ~ exp(— ip v x v ) where p v is the 4- 
momentum. Imposing the Lorentz guage condition d^U^ = demands p^e^ = 
0. Therefore the number of independent polarization vectors is reduced from 
four to three in a covariant fashion. However, one can still perform another 
gauge transformation to the massless U® which makes finally only two degrees 
of freedom remain. Therefore, one should keep in mind that in the present 
model the spatial velocity has only two degrees of freedom, that means one 
component must be described by another two vector components. Fortunately, 
in real applications in cosmology or compact star, this requirement is satisfied 
by the assumption that the system under consideration is isotropic. 

On the other hand, let us comment on the gauge transformation of velocity 
m£. The gauge field [7° is transformed as U% -> (U£)' = U% + l/g s df,0 a + 
jmbcQbjjc w j^h the local gauge parameter 9 a = 6 a (x). Following the Lorentz 
condition, 9 a also satisfies d 2 6 a = 0. If the solution is 9 a = in a (f)cxp(—ip ■ x) 
with n a constant and <p ~ cxp(— • x), is transformed as u 1 ^ — » = 
u^ l + K a /g s p t t+if abc K b u'! 1 . The invariant 4- velocity satisfies u' 2 = l+n(m g / 'g s + 
if g ) + 0(k 2 ) ~ 1 = u since the second and third terms are suppressed by 
K a ~ 9 a <C g s ~ 0(1). f g is the factor of summed colored gluon states. 
Therefore, one can take the form of Eq. (2) for a good approximation. 

From now, throughout the paper let us focus only on the gluon sea of 
plasma. This means one should consider only the related gluonic terms in Eq. 



3 Energy momentum tensor 

Now we are ready to proceed with deriving the energy momentum tensor 
within the model. It should be pointed out that once the hot (high energy) 
QGP state is achieved, the system is assumed to be predominated by the 
classical motion rather than the quantum effects. 

Therefore the total action of matter for non-gravitational fields in a gen- 
eral geometry of space-time 1Z is S g = J n d 4 x ^/^g£ g , where g is the deter- 
minant of metric g^ v . The variation of S g in the metric is given by SS g = 
~\ In ^ x V^T^v 5g^ u . Since the energy momentum tensor density is, 



(1) 




na na^v . jajja^ 



(5) 




2 8C g 



(6) 
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one obtains, 

= S® p S p v — g^ v Cg + 2g s J^U v . (7) 

It is clear that Eq. (7) is symmetric as expected to fulfill the Einstein gravita- 
tional EOM. The total energy momentum tensor T M „ is given by integrating 
out Eq. (7) in term of total volume in the space-time under consideration. This 
means T^ v is a result of bulk of gluons flow in the system. 

Furthermore, in a general space-time coordinates, the components of en- 
ergy momentum tensor determine the total energy density (Too), the heat 
conduction (To^jo), the isotropic pressure {To) and the viscous stresses (Ty 
with i 7^ j) of the gluonic plasma. Of course, in this case the derivative in- 
side the strength tensor should be replaced by the covariant one, V„. Also, 
the energy momentum tensor satisfies the conservation condition, V M l~^ v = 0. 
Nevertheless, one can trivially conclude that the model induces non-zero vis- 
cosity since generally =/= for i =/= j. From the experimental clues, however 
the size should be small such that it is always treated perturbatively in most 
hydrodynamics models [8-13]. 

Before going further to apply these results, one should determine the quark 
current J° in Eq. (7). This can be simply calculated by considering the EOM 
(Dirac equation) of a single colored quark (Q) or anti-quark (Q) with 4- 
momcntum p M . Since the solution of the EOM is Q(p, x) = q{p) exp(— ip ■ x), 
one immediately gets q^^q = 4p M . Assuming that all colored quarks / anti- 
quarks have the same momenta and the velocity of gluons and quarks inside 
the gluon sea are homogeneity, approximately J®U a ' J ' oc Ap^U 11 = AmQ<j) since 
= u 2 = 1. 

4 Equation of state 

Now we are ready to consider the compact star interiors in the model, partic- 
ularly before transforming itself into neutron star. 

The whole compact star is commonly described as a static spherically sym- 
metric space-time represented by Schwarzschild geometry. This means one 
deals with the relativistic gravitational equations for the interior of spherically 
symmetric plasma distribution in the core. In the region under consideration 
the presence of gluonic fields flow induces non-zero energy momentum ten- 
sor which is making up the star. This is the phase before the neutron star is 
getting mature. Starting from the stellar nebula made of hot plasma which is 
gradually getting colder as the hadronization occurs from the colder surface, 
while the inner core is still in pure hot plasma state. 

As a consequence of the diagonal metric of Schwarzschild space-time, the 
model falls back to the perfect fluid without viscosity and heat conduction, 
i.e. T 0i = Tij = for i ^ j. Also, since the plasma distribution should 
be spherically isotropic, it is considerable to put v\ = V2 = t>3 = v as con- 
stant for all colored gluons. This assumption is consistent with the degree of 
freedom counting discussed in the preceding section. Moreover, the vanishing 
off-diagonal components of the Ricci tensor, R i0 , actually forces the spatial 
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3-velocity of the fluid must vanish everywhere. Hence particular assumption 
for Vi is indeed not necessary. However, the gluon distribution still depends on 
the radius length, = 4>{r). 

For the sake of simplicity one can put homogeneous gluon fields for all color 
states, i.e. = for all a = 1, • • • , 8. This yields, 



% v = [8 g s f Q m Q <f>(r) + g 2 s f g <f>{r) 4 ] u^u„ 



4,g s .f Q m Q ^r)--g 2 f 2 0{r) 4 



9tiv ■ 



(8) 



where f g is the factor of summed colored gluon states from the structure 
constant f abc , while /q is the factor of summed colored quark states from 
jajjan R emm( j that the energy momentum tensor for perfect fluid takes the 
form, 

T flv = (£ + V)u fi u v -Vg^. (9) 

Here £ and V denote the density and isotropic pressure for single fluid field, 
each is related to the total density and pressure of the system through p = 
§ d 4 x£ and P — § d 4 xV respectively. Obviously, from Eqs. (8) and (9) one 
can obtain the density and pressure in the model as follows, 



P(r) 



p{r) 



f" dt [ dV 

= ^9sfQm Q 
T 

= f' dt f d 



1 



45s /Q m Q (f>(r) - -g s f 4>{r) 



dV 



9s ft 



16 /q toq 



<t>{rf 



(10) 



^9sfQm Q 



/ 



dV 



1 + 



16 /q m Q 



4>{r). 



(11) 



at a finite temperature /3 = 1/T in a 3-dimensional spatial volume V. Here, 
T s and To denote the core surface and inner core temperatures. 

The proper spatial volume element for Schwarzschild geometry is dV = 
\J B(r)r 2 sin 9 dr d9 dip with radius r and two angles 9 and <p in spherical 
coordinates. The solution for B(r) is given by B(r) = [1 - 2Gm(r)/r] with 
m(r) = 4ir JJ dfp(f) f 2 is the 'bare mass'. This generates the proper integrated 
mass m(r) contained within a coordinate radius r inside the star. Anyway, the 
inner structure of star with Schwarzschild geometry is well known as Tolman- 
Oppcnhcimcr-Volkoff (TOV) equation which relates density and pressure in a 
unique way [32,33]. 

Obviously, using Eqs. (11) and (10) one can construct certain EoS, 



P(r) = w{r) p(r) , 



(12) 
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where in contrast with the conventional cosmological models, 



w(r) 



I 



dV 



I 



dV 



f 2 



16 Jq rriQ 



0(r) 3 



1 + 



16i^ 0(r) 



<f>(r) 



(13) 



<f>(r) 



is not a constant. This is actually one of the important consequences in the 
present model. Once the field distribution <\>(r) is determined one can obtain 
certain forms of density and pressure. 



5 Pressure and density distribution 

For the sake of convenience later on, let us redefine the radius r to be the 
dimensionless one, i.e. the ratio of core and compact star radius : r — » r' = 
r/r . Here r is the compact star radius. Then, the pressure and density in 
Eqs. (11) and (10) can be expressed as, 



F x (r ) = F Xo + 16tt g s f Q m Q r a 



( dr'r' 2 ^B{r') 
Jo 



3 T — T s 
T T S 

l + k x 



9sf 9 . (/)3 



16 fQ m Q 



4>{r') , (14) 



after integrating out the time component. Fx denotes the pressure P or density 
p for X = P, p respectively, while kp — — 1 and k p = 5. Fxo represents the ini- 
tial F x . Taking A^r') = J^' dr'r' 2 ^/B{?)<f>(r'), A 2 = ^ dr'r' 2 /Bp)<£(r') 4 , 
ki = 16tt g s f Q m Q r a 3 and k 2 = tt g 2 s f g r Q 3 yield, 



T n -T s 



Fx(r') = Fx + ^ F 7r L [kiA 1 (r') + kxk 2 A 2 (r')] . (15) 



Finally, by substituting the expressions for m(r') and B{r') into Eq. (15) 
one finds, 



■p / i\ el T Q - T s r' [k 1 + kxk^jrf] 
F x (r)=F X0 + / dr = 

JO-'s Jo 



(16) 



f dr'p(r')r' 2 
r Jo 



Taking its derivative in term of r' , Eq. (16) leads to an ODE as below, 



^(rO^^O+r 2 ^^) 3 -^ 2 ^^^)^^) -^^')^^) = , (17) 



-2/2 
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where T = T T S /(T - T s ) and, 



A x (r>) = 2r'^(r') 2 [h + k x k 2 ^r') 3 } 2 , (18) 
A 2 = 8nG , (19) 

A 3 (r') = [$r'U(r r f + 2r'^{r')4>' (r 1 )] [h + k x k 2 ^(r'f] 2 

+6k x k 2 r' 5 [fci + k x k 2 4>{r'Y] 0(r')V (r') . (20) 

Note that a prime in <j> or F x means a derivative in term of r'. 

Nevertheless, one can also derive a "differential" EoS by taking the deriva- 
tive of Eq. (16) in term of r' and eliminating B(r'), that is, 



l(>) = l&f Q m Q -g s p g cp{r'f 
K > 16/ Q m Q + 5 3s /^(r') 3Pl >' 



(21) 



Now let us solve the ODE in Eq. (17). The ODE can be solved analytically 
by taking an approximation of small compact star core, i.e. r'<l. This ap- 
proximation is quite natural since by definition the compact star cores should 
be small enough. Hence one can expand the ODE near the origin (r' — > 0) to 
obtain the solution order by order, 



F x (r') ~ F x (r>) + F^ (r')r' + ^F^ (r')r' 2 + ±4"V)r' 3 + ' ' ' 

First of all, taking the first derivative on Eq. (17) in term of r' yields, 

3T 2 F x I)2 (0)F x n) {0) = , 
since A'^O) = A' 3 (0) = 0. Further derivatives up to 7th order give, 

1 \2A 2 F x (0)FP(0f - 64 7) (0)4 77) (0) 21 



F X IH \0) 



3F x I] (oy 



r'=0 

(22) 



(23) 



(24) 



F X IV \0) 



F X V \0) 



1 



F { x \0f 



6A 2 F X (0)F X I) (0) 2 F X II \0) + 2A 2 F X 1 >(0) 



(')/ 



-2T 2 F ( x I] (Qf - 6F^ 7) (0)i^ 77) (0)F^ 777) (0) 



(25) 



3T 2 F^ 7) (0) 2 



,, \Q)F^ ) {{)) + 'dQA 2 T 2 F x {Q)F' x 1 '{QfF^ ll '{Q) 



+72yl 2 T 2 F x (0)i^ 7) (0)F^ 77) (0) 2 + 82/l 2 T 2 F^ 7) (0) 3 i4 77) (0) 
-2AT 2 F ( x \q)F ( x I] {{))F { x V \q) - 18T 2 F^ 7) (0)i^ 777) (0) 2 
-36T 2 ^ 77) (0) 2 ^ 777) (0) 



(26) 



4 y/) (o) = 



3T 2 F^ 7) (0) 2 



7&A 2 T 2 F ( x \q) 6 F 



3^(111), 



x (0) + 36A 2 T 2 F x (0)F ( x I) (0f 



+24A 2 T 2 F X (0)F ( X \0) 2 F ( X IV) (0) + h22A 2 T 2 F { x \{)) 2 F x II] {Q) 



0.2 0.4 0.6 0.8 1 




Fig. 1 The density distribution as a function of the normalized compact star core radius 
with T a = 175 MeV and T = 1 GeV up to 4 th (dotted-dashcd line), 5 th (dashed line), 6 th 
(thin solid line) and 7 th (thick solid line) order expansions. 



+84/l 2 T 2 F x (0)i4 //) (0) 
+M2A 2 T 2 F X (0)F { X I] (0)F X II) 

(IV 



? {VII) 



(0) 



188,4 2 T 2 fi 7) (0) 3 F (777) 



(o)4 7// 



+36A 2 T 2 F X (0)F^ 7) (O) 2 ^ v ' (0) + A ( 3 V > (0)F X J > (0) 



x (or 

>(0) 



(0) 



+ 5Ai IV \0)F^(0) 



(V) 



4-) ( o) 



72T 2 F X II> (0)F X ' 11) (Q) 



-60T^ 77) (0) 2 F^(0) - 18T^ JJ (0)Fi JJJj (0) 



(ii), 



x 
pi 111 ), 



(in), 



-60T 2 F^ 7) (0)F^ 77) (0)F 



(IV), 



x (0) 



-30T 2 F { X ] (0)F x n) (0)i4 y) (0) 



3T 2 F { x \oy 



x 

378A 2 T 2 F X 



(27) 



2 ^ (7) (0) 3 4 7y) (0) 



4 W) (0)F 



+66 A 2 T 2 F X (0)4 7) {0) 2 F X V) (0) 
+154 7y) (0)4 77/ )(0) - 6A^(0)F X IH \0) 



(i) 



(0) 



-72T^4 77) (0) 4 - 18T^ JJJj (0) 3 - 60™^(0)F^(0) 



,(/// 



(i) 



?(iv). 



-90T 2 4 7) (0)4 777) (0) F x } (°) 



(28) 



From Eq. (23), non-trivial solution is obtained for 4 77 ^(0) = which 
leads to F x \o) = const. Without loss of generality, one can choose Fx(0) = 



60 




0.175 0.225 0.275 0.325 0.375 

T (GcW) 

Fig. 2 The density distribution as a function of inner temperature with the normalized star 
radius r' = 0.01 and T s = 175 MeV. 



1. These results are substituted subsequently into Eqs. (24)~(28) to get the 
complete solution. The solution is depicted in Figs. 1 and 2 for the density 
distribution, i.e. Fx = p. 

It is important to derive the solution analytically to investigate the be- 
havior of phase transition in term of higher order solutions. The results of 
total density as a function of r' arc depicted in Fig. 1 for the solution up to 
4 th (dotted-dashed line), 5 th (dashed line), 6 th (thin solid line) and 7 th (thick 
solid line) order expansions. It is clear that the phase transtion occurs at 7 th 
order accuracy. In particular the contour is mainly governed by the dissipa- 
tive term, A^ T \o). Those results are obtained by taking a particular form of 
distribution function, <p(r) ~ exp(ar) with a = 0.4. 



6 Summary 

The QGP, in particular gluon matter, dominated compact start core interior 
has been investigated using the fluid QCD lagrangian. From the lagrangian 
one can derive the energy-momentum tensor and subsequently the density and 
pressure distributions from the first principle. In the model, both density and 
pressure are related each other and form an extended EoS of perfect fluid 
through the field distribution 4>{r). The relation leads to the linear EoS known 
in the conventional cosmological model, but allowing the EoS parameter w to 
be a function of temperature T and core radius r. The analytical solutions for 
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the total pressure and density distribution equations are also given for small 
normalized compact star core radius, r'<l. 

The model describes the dense gluon plasma inside the compact star core as 
a bulk dynamics of gluonic matters before the cooling phase near the compact 
star core surface as depicted in Fig. 1. It should be noted that the figure 
provides a general contour for pressure or density distribution inside the core. 
This suggests that near its surface the hadronization is getting involved to 
form surfaces of small compact star cores, 

On the other hand, the phase transition for certain size of compact star 
core in term of temperature changing is shown in Fig. 2. It can be observed 
that the phase transition occurs at the scale of hadronization around the core 
surface, while it keeps fluidity as the temperature is getting higher in the core 
interiors. 
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